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Abstract: We show that the thermal rehc abundance of dark matter can be affected by a new 
type of reaction: semi-annihilation. Semi-annihilation takes the schematic form tpitpj — )• ipk^P^ 
where ipi are stable dark matter particles and is an unstable state. Such reactions are 
generically present when dark matter is composed of more than one species with "flavor" 
and/or "baryon" symmetries. We give a complete set of coupled Boltzmann equations in 
the presence of semi-annihilations, and study two toy models featuring this process. Semi- 
annihilation leads to non-trivial dark matter dynamics in the early universe, often dominating 
over ordinary annihilation in determining the relic abundance. This process also has impor- 
tant implications for indirect detection experiments, by enriching the final state spectrum 
from dark matter (semi-)annihilation in the Milky Way. 
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1. Introduction 

The existence of dark matter is one of the best motivations for physics beyond the standard 
model (SM). Evidence for dark matter has accumulated at vastly different length scales — from 
galactic scales and clusters of galaxies to global scales of hundreds of megaparsecs [Q, ^ ^ . 
However, all of these observations infer the existence of dark matter through its gravitational 
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effects alone. In particular, they do not tell us the nature, origin, or composition of this 
important component of our universe, which is not explained by any SM degree of freedom. 

A particularly well-motivated class of dark matter candidates are so-called Weakly Inter- 
acting Massive Particles, or WIMPs, whose abundance is determined through thermal freeze- 
out. In the Lee- Weinberg scenario |^], WIMPs ifj thermalize in the early universe through 
the annihilation reaction ipip — )• where (/> and (/>' are SM degrees of freedom, until their 
interactions freeze out. A standard relic abundance calculation ||5|, ^, ^, ^, |9|, |ll| shows that 
the dark matter mass density today depends only logarithmically on the relic particle mass 
and scales inversely with the total annihilation cross section at freeze-out, QnMh'^ oc {av)^^. 
It is remarkable that for a dark matter mass between 10 GeV — 10 TeV and an electroweak 
annihilation cross section, one gets a relic abundance in the ballpark to account for dark 
matter. This fact is often referred to as the "WIMP miracle" , and in the context of the gauge 
hierarchy problem, it is suggestive that the same particles one might introduce to stabilize 
the Fermi scale could also explain the dark matter in the universe. 

Since the WIMP paradigm is one of the best explanations for dark matter and since 
thermal freeze-out is so predictive, it is important to know how to correctly compute the 
dark matter thermal relic density. This is particularly important for models where the relic 
computation cannot be reduced to the Lee- Weinberg scenario. As we will argue, the thermal 
abundance of dark matter can be dramatically affected by the presence of a new dark matter 
interaction, which we call "semi-annihilation". 

Semi-annihilation occurs when dark matter is stabilized by a larger symmetry than just 
Z2. In the simplest case with just one dark matter species ip, there can be an additional 
allowed reaction 

^p^p — )• tp(p (1-1) 

which preserves a symmetry. Here, i;^ is a SM state or a new particle which decays to the 
SM. We see that unlike ordinary annihilation where the total dark matter number changes 
by two units, in semi-annihilation the total dark matter number changes by only one unit. 
More generally, dark matter can be composed of more than one stable component ipi, and 
these relic particles can have non-trivial mutual interactions. In this case, a more general 
semi-annihilation reaction is possible, 

ipiipj tpkcl), (1.2) 

which often occurs if dark matter is stabilized by "baryon" and/or "flavor" symmetries, as in 
QCD-like theories. Such reactions are also allowed in models where dark matter is composed 



of non-Abelian gauge bosons ||l2|, 13, In this paper, we study explicit examples of such 
models and find that to correctly compute the relic abundance, semi-annihilation must be 
included. In fact, semi-annihilation can dominate over standard annihilation for some regions 
of parameter space. 

To understand how semi-annihilation fits into the WIMP paradigm, consider a more 
general framework for dark matter interactions. Start by assuming the existence of a new 
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Table 1: Different dark matter freeze-out scenarios and allowed reactions. We do not include the 
reaction ipitpj — > ipk'4'm in this table since it is always present, and we explicitly indicate when only 
the diagonal (z = j) contribution is allowed. 



dark sector composed of N particles ■0* which can be either stable or unstable. The possible 
reactions involving ipi which can take place in the early universe are 

a) tpiipj <j)(j)', b) Tpi4> ipj<j)', c) ipiipj Tpkcj), d) ipiipj tpkipm, e) tpi ipjcj), (1.3) 

where again (j) and cp' are light degrees of freedom in thermal equilibrium with the SM. 
Different dark matter thermal freeze-out scenarios depend on which of the reactions a)~e) are 
active, and we summarize these main possibilities in Table 

Lee- Weinberg: The simplest case is the Lee- Weinberg scenario [^, where = 1 and the 
only allowed reactions are a), b), and d), with i = j. Chemical freeze-out is determined 
by the annihilation reaction a), and kinetic freeze-out is determined by b). Reaction d) 
plays no role in the thermal relic computation. 

Co-annihilation: A slight variation of the standard case is when > 1 but there is 
only one stable dark matter species. Here the relevant reactions are a), b), d), and e). 
This is the case in the minimal supersymmetric standard model with i?-parity, where 
all heavier particles are unstable and eventually decay to the lightest one via reaction 
e). In principle, this co-annihilation case involves a system of coupled Boltzmann 
equations, but as we will review, as long as reactions of type b) with i ^ j are effective 
at freeze-out, it is possible to compute the relic density via standard methods llf 



Decoupled Mult i- Component: These first two examples assume dark matter to be com- 
posed of a single particle, but more generally, dark matter could be composed of more 
than one stable component. Many such multi-component dark matter models have 
been proposed (see e.g. |16, 17, 18, 18, 19, |0|, |2^). The standard approach in 



multi-component models is to assume that each particle thermalizes independently of 
the others, thus the total dark matter density today is JIdm = where the sum 

runs over all the thermal relics. Said in the language of Eq. ( |1.3D , only the diagonal 
(i = j) reactions of type a) and b) are present, and reactions c) and e) are forbidden. 
Reaction d) may or may not be present in such models, and is relevant for calculating 



the relic density [21|. 
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Semi-annihilation: The focus of this paper is on the reaction of type c) , which to our knowl- 
edge first appeared in Ref. |12, |l3|, 0].^ As long as the triangle inequality rrik < mi+nij 
is satisfied (as well as its crossed versions), the semi-annihilation reaction V'iV'i ~^ i^k4' 
can take place without making any relic particle unstable. Semi-annihilation implies 
that the relic particles have non-trivial mutual interactions, and typically both reac- 
tions of type c) and d) are important in determining the dark matter relic abundance. 
For simplicity, we focus on the case where each particle in the dark sector is absolutely 
stable, in which case reaction e) is forbidden. In addition, only the diagonal reactions 
of type a) and b) are allowed, since off-diagonal contributions would imply that the 
heavier particle would be unstable by crossing symmetry. As we will see, the absence 
of off-diagonal reactions of type b) makes these models more difficult to study than the 
standard cases. 



The remainder of this paper is structured as follows. In Sec. § we present a toy model 
where dark matter is composed of a single component but semi-annihilations are present 
and can dominate dark matter production in the early universe. In Sec. |^ we introduce 
example multi-component models where semi- annihilation reactions are present, and we give 
a complete set of Boltzmann equations which can be solved to obtain the relic density today. In 
particular, we argue that a semi-analytical solution analogous to the Lee- Weinberg scenario or 
co-annihilation is not in general possible, and thus the equations have to be solved numerically. 
We present a minimal multi-component toy model with semi-annihilations in Sec. Q and 
numerically study the effects of semi-annihilation on the relic abundance. We explore the 
effects of semi-annihilations on indirect detection experiments in Sec. ^, and conclude in 
Sec. p. Computational details, as well as an example supersymmetric QCD model with 
A'^^ = Nc + 1, are given in the appendices. 



2. Semi-annihilation with a symmetry 

To introduce semi-annihilation, we present a simple case where semi-annihilations play a 
significant role in the thermal relic calculation. We assume dark matter to be composed of 
one single stable particle ^ complex scalar with mass m^, which is stabilized by a Z3 
symmetry.^ The x particle interacts with a real scalar (/), which eventually decays to SM 
states. The interaction Lagrangian is 

= ml x^X + ai X^X(t> + a2 X^X<P'^ + ^3 - ^(0), (2-1) 



^In the models considered in Ref. |12| |l^, [lij, a custodial symmetry effectively reduced the relic abundance 
computation to a single particle system. Here, we consider semi-annihilation in multi-component dark matter 
models where there is no such simplification. 

^For a concrete realization of dark matter models with a Z3 symmetry and with semi-annihilation playing 
a subdominant role see Ref. |24[ H]. See also Ref. jie) and Ref. |||, H. 
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Table 2: Field content and symmetries of the model with a Z3 symmetry. 
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(a) 



(b) 



Figure 1: The reactions that keep x in thermal equilibrium in the model with Z3 symmetry: (a) 
annihilation; (b) semi-annihilation. The x fields are drawn as double lines. The field decays to SM 
states such that the whole system is in thermal contact with the SM before freeze-out. 



where contains additional interactions for ^ alone. This function contains both a scalar 
potential for ^ and couplings with SM fields, 

V{(\)) = ci (/.^ + C2 </>^ + di (PBf.^B"" + d2 (t>H^H + • • • , (2.2) 

where B^y is the hypercharge field strength, H is the Higgs doublet and the • • • stand for 
other possible interactions. The model and its symmetries are summarized in Table ^. 

As long as the coupling strength of (j) to SM fields is sufficiently large, (p will be kept in 
thermal equilibrium with the SM. Then, the x particles are kept in thermal equilibrium in 
the early universe by the reactions XX ~^ (t"P ^-nd XX ~^ X0i ^-s shown in Fig. |l|. While we 



focus on the case of a field <p "portal" to the SM |2£, 31, 32], in more general dark matter 



scenarios (p could be a SM field itself |^3|, |3j, |35| . 
2.1 Boltzmann equation 

The Boltzmann equation describing the number density evolution of x is^ 



- + 3-?^% = -(^^)xx^00 [4 - "^7^] - (^^)xx^x</' hx - "^x^^] ' (2-3) 



dt 



■^For a brief derivation of the Boltzmann equation in an expanding universe, see App. |^. The following 
Boltzmann equation also appears in Ref. [ |l3| , , despite the fact that they have a multi-component model. In 
that case, a custodial symmetry implies a degeneracy for the dark matter species number densities, reducing 
the relic abundance calculation to a single differential equation. 
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where H is the Hubble parameter and n'^ is the equiUbrium number density distribution. 
In the (crv)^^_j.i^0 = hmit we recover the famihar Lee- Weinberg scenario. To get an un- 
derstanding for the behavior of this Boltzmann equation, we restrict ourselves to s-wave 
(semi-)annihilation, and we choose as free parameters the s-wave amplitudes for the two 
processes, rja and r]s respectively. The thermal averages of the s-wave cross sections can be 
obtained using Eq. ( |C.5| ) from App. |C|:^ 

A, A, 

2.2 Semi-analytical solution 

The above Boltzmann equation can be solved semi-analytically in analogy with the Lee- 



Weinberg calculation. We rewrite Eq. ( |2.3| ) by introducing a dimensionless time variable 
X = m^/T and comoving number density = n^/s, where s is the entropy density of the 
relativistic degrees of freedom. The Boltzmann equation in these new variables is 

^=4k-';'v^k-w]. ,2.5, 

We can give a semi-analytical solution to this equation following the method in Ref. by 
solving the equation in two regimes, early and late times, and then matching the two solutions 
at the freeze-out point. For convenience, we define the function A = 5^ — Y^. 

At early times, the number density is very well approximated by its equilibrium value, 
thus we impose dA/dx = in the Boltzmann equation. The use of the Maxwell-Boltzmann 
equilibrium distribution for Y^"^ is well-justified for a cold relic. Eventually a point de- 
noted "freeze-out" is reached when the interaction rate is not fast enough compared with 
the expansion rate to maintain thermal equilibrium. The freeze-out point Xf is defined as 
A(xy) = cY^'^{xf), where c is a numerical constant of the order 1. Since we are making 
the s-wave approximation, we take c = \/2 — 1 as suggested by Ref. After freeze-out, 

the equilibrium distribution is exponentially suppressed and we can safely neglect it in the 
Boltzmann equation and solve directly for Y^. We finally match the two solutions at the 
freeze-out point Xf. 

The freeze-out value Xf is found by solving the equation 



Xf = log 



0.038 c(c + 2){av)^^^^^^^^-^J^ 



+ log 



c + 2 (av) 



XX- 



(2.6) 



where = 1 for a scalar field, g^, is the effective number of relativistic degrees of freedom at 
the time of the freeze-out, and the Planck scale is Mpi = 1.22 x 10^^ GeV. The first term on 
the right-hand side of Eq. ( |2.6| ) is the solution we usually have in the Lee- Weinberg scenario. 
Thus we see that the effect of semi-annihilation is to shift the freeze-out temperature by only 
a small logarithmic amount. 



*The factor of 3/4 comes from the phase space suppression in semi-annihilation. 
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The mass density of the x particle today is given by p^^o = soY^{oo), and as usual 
can be expressed as the fraction of the critical energy density 

0,.'=2xl^Iiii5!^, J(.,)= rrf.(^!;>M^i»±MM^, (2.7) 

where the 2 factor takes into account the antiparticle degree of freedom and the annihilation 
integral includes both annihilations and semi-annihilations. 

2.3 Results 

The results of solving Eq. ( |2.5| ) are shown in Fig. || for the mass value = 1 TeV. We 
consider three representative cases: annihilation only, semi-annihilation only, and both pro- 
cesses present with equal amplitudes. In each case the parameters are chosen in order to 
get the observed relic abundance, VL^h? ~ 0.1 [^6|. We plot both the numerical and the 
semi-analytical solutions, and we find a good agreement between them. In Fig. §(d), we show 
the ratios between the various cases, where we see that as we increase the semi-annihilation 
contribution, it takes longer for the system to reach the freeze-out density. This behavior 
is expected from Eq. (|2.5| ), since the collision operator for semi-annihilation has a positive, 
linear dependence on that somewhat counteracts the overall depletion of x particles. 

In Fig. ^, we identify the region of the {r]a,ris) plane which gives the observed dark 
matter relic density in our universe. We consider two different dark matter particle masses, 
rriy. = 1 TeV and niy. = 5 TeV respectively, and we shade the region where the dark matter 
relic density is within the WMAP 95% CL region namely 0.0975 < ^dm/i^ < 0.1223. 
In single component dark matter models, the conventional assumption is that we have only 
annihilation (namely 7/^ = 0), thus once we fix the dark matter mass niy. the amplitude r]a 
is uniquely determined by the relic density measurement. However, if we allow also semi- 
annihilations, a wider region of parameter space gives the correct relic density. Such a region 
corresponds roughly to At]'^ + 3r]^ ~ 1.5, thus the thermal production of the relic particles in 
the early universe can also be completely controlled by semi-annihilations. 

3. Semi-annihilation with multiple species 

The previous toy model shows that semi-annihilation does affect the dark matter relic density. 
However, single-component dark matter with a symmetry is not representative of the 
kinds of models in which semi-annihilation is relevant. In this section, we introduce an 
example multi-component model having semi-annihilations among its allowed reactions. A 
more realistic construction is given in App. ^ where a supersymmetric gauge theory with 
Nf = A^c+1 is considered. Here, we prefer to focus on a simpler model with a similar symmetry 
structure, thus avoiding unnecessary complications caused by having superpartners. 

3.1 A meson-baryon system 

The dark matter model we consider consists of fermions and bosons that have a "baryon" 
number symmetry as well as a large "flavor" symmetry. The system contains A'^^ dark baryons. 
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i]a=0.61, i]s=0, 0.h^ =i. 0.1 



j/„=0, /7s=0.72, 0.h^ 0.1 




(c) (d) 

Figure 2: Comoving number density Y{x) in the x niodel witli a Z3 symmetry for my- = 1 TeV. Three 
cases are considered: (a) pure annihilation, (b) pure semi-annihilation, (c) equal amplitudes for the 
two processes. In all the three cases, the parameters rja and r/s are chosen in order to get the WMAP 
value for the relic density today, flY^uh^ ~ 0.1. We plot the equilibrium solution (dashed black), the 
numerical solution (solid blue) and the semi-analytical solution (dashed green). To emphasize the 
subtle differences between the cases, in (d) we show the ratio between each numerical solution and the 
one for pure annihilation, where pure annihilation is yellow, pure semi-annihilation is blue, and the 
mixed case is green. 



namely vector-like fermions 6* and 6?, where z is a flavor index. ^ In addition, there are Nj 
dark mesons, namely complex scalars xl arranged in an adjoint of U{Nf). The charged dark 



baryons and mesons interact with a neutral portal field = Tr 



as in Sec. 0. The particle 



content and symmetry structure is summarized in Table H. The Lagrangian of this system is 



-bbx 



m,,^b^b' + Xb''^xlb' + h.c. -Vixi), 



(3.1) 



^We work in two-component notation, where b and b'^ are independent left-handed Weyl spinors. 
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m^ = l TeV m^ = 5 TeV 
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(a) (b) 

Figure 3: Relic density for tlie x model with Z3 symmetry in the {r]a,T]s) plane for relic particle 
mass: (a) = 1 TeV; (b) = 5 TeV. The shaded regions correspond to 0.0975 < fl^h'^ < 0.1223 
(WMAP 95% CL region [Q). The dashed lines corresponds to constant dark matter density, and 
they are along the curves A-q^ + irjl ~ 1.5. 
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Table 3: Field content and symmetries of the meson-baryon system. 



where the scalar potential V{Xi) contains all the renormalizable operators involving the meson 
matrix field consistent with the flavor symmetry in Table ^. In particular it contains the 
interactions for the singlet (j) = Tr xi analogous to the ones in Eq. (p.2|), allowing such a 
singlet to decay to SM particles. 

In the limit of an exact U{Nf) flavor symmetry, all the particles except (p are charged, 
therefore all of the mesons are degenerate and all of the baryons are degenerate. In addition, 
if the triangle inequality rule is satisfied for the decay xl ~^ hbj, then the charged mesons and 
baryons are stable, and we can consider all of them as dark matter candidates. If a spurion 
\l is introduced, with (A;-) = SijXi, the original U{Nf) flavor symmetry is explicitly broken 
to [/(l)^^. The lightest baryon is still stable by the U{1)b symmetry. In fact, the heavier 
baryons might be stable as well. While the flavor symmetry allows 6* to decay to the lighter 
state bj through the process ¥ — )• bjxl ! such a decay is forbidden by kinematics if the triangle 
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Figure 4: Triangle inequality rule for the stability of particles tpi, tpj, and ipk- In the case on the left, 
the masses can be put on a triangle and thus the three particles are mutually stable, whereas for the 
case on the right, this is not possible and the decay ipk — >■ i^i4'j4' is kinematically allowed. 



inequality rule mi,. < nib + m j is satisfied (see Fig. 0). The off-diagonal mesons {i ^ j) 
are still charged, thus the lightest one is absolutely stable. The heavier off-diagonal mesons 
might be stable as well if their symmetry-allowed decays are kinematically forbidden, namely 
if they satisfy the triangle inequality. Note that the diagonal mesons do not carry any flavor 
quantum numbers, thus there is no symmetry which guarantees their stability. We denote 
them by xl = and the way they decay to SM fields is through operators analogous to 
Eq. (P). 

It is now clear from the symmetries of the theory that the semi-annihilations reactions are 
indeed present when the stable dark matter particles are thermally produced at the freeze-out, 
as for example 

x\xl^xl<P\. h^W^hi<p\. (3.2) 

We now give a complete set of Boltzmann equations which can be solved to get the relic 
abundance in the case of semi-annihilations, and we also explain why in general a semi- 
analytical solution is not possible in this case. 

3.2 Boltzmann equations 

We consider a generic dark matter model with stable components and assume that the 
symmetry structure of the theory allows semi-annihilations to take place at the freeze-out. 
We ignore unstable states in the dark sector, though in principle they would change the 
relic density through co-annihilation. Given our assumptions, the reactions which have to be 
included in the Boltzmann equations are the ones in the last row of the Table |l|. If we label 
the stable particles as V'ij where the index i goes from 1 to iV, the Boltzmann equation for 
the number density rii of the tpi particles result in 

^ + ?,Hni = - {auVj-ci) (nf - n^"*^^ - '^{aijkVrci) (^iUj - ^n^n^ 



eq eq ^ 



(3.3) 



j,k,m \ 



m 



where an is an ordinary annihilation cross section ipiipi — )• (jxj)' , fijjfc is a semi-annihilation 
cross section ipiipj — )• V'fc^j and (Jijkm is a dark matter conversion process ^pitpj — )■ tpki^m- The 
collision operators on the right-hand side are derived in App. |A| . 
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The pure annihilation contribution an must be diagonal in dark matter flavor space, 
since if it were not, the heavier particle between ■0j and ipj would be unstable by crossing 
symmetry. As long as the triangle inequality rule from Fig. ^ is satisfied, then the semi- 
annihilation process aijk has no such restriction. The reason why we need to include dark 
matter conversion (Jijkm is because we care about the total energy density p = niUii, and 
since all of the particles are stable, such reactions can change p if the particle masses are 
different. Finally, no reaction of the type V'i'/' — ^ ipj'P' takes place, except for the diagonal one 
i = j, by the same stability argument, and we do not have to take into account the diagonal 
process since it does not change the total number of ipi particles. The lack of the off-diagonal 
contributions for this process is the reason why in general we cannot have a semi-analytical 
solution to the system in Eq. ( |3.3D , as we now show. 



3.3 Comparision with co-annihilation 

The system in Eq. (3.3) contains coupled equations of the Riccati type, thus it cannot 



be solved analytically. In the case = 1 it is possible to solve the single equation in two 
different regimes, match the solutions at freeze-out and get a semi-analytical result for the 
relic density today, as done in Sec. |2|. 

However it is well known that in the co-annihilation case, described in the second row 
in Table ||, we can sum the equations in Eq. ( |3.3[) to get a single one of the same form as 
the Lee- Weinberg scenario, even if at the beginning we were dealing with a system of A^ 
coupled Boltzmann equations [^]. The crucial difference between co-annihilation and semi- 
annihilation scenarios is the presence of off-diagonal ^picp — )• V'i^' reactions in the former. Such 
reactions are much more effective at freeze-out than V'iV'j (co-)annihilation, since to have co- 
annihilation, we need two dark matter particles in the initial state and their number density 
is Boltzmann suppressed at that time, whereas the <j) particles have a relativistic number 
density. 

A more quantitative way to see this fact is to write down the collision operator for the 
—7- ^j</>' reaction 



eq 

n- 



C\l}i<l>^i)j(t)' — nin^{av)^.^^^^^i 1 |q . (3.4) 

V J 

Since there are a lot of ^ particles around, this reaction is very effective and it guarantees the 
condition 

^ = !C, = (3.5) 

rij Tij n n^^ 

Thus in the case of co-annihilation (but not in semi-annihilation) we can assume the condition 



in Eq. (3.5) before, during, and after freeze-out. This relation between equilibrium and 
non-equilibrium number densities allows for a considerable simplification of the Boltzmann 
system. Summing these equations together, we get a single equation governing the total 
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number density n = Y^- ni involving an effective annihilation cross section a^s 1 15] 



— + 3Fn = -{UeSVrel) {n^ - rf^"^) , (Jeff = ^ <Tij ViVj, (3.6) 

i,3 

where aij is the cross section for the V'iV'j co-annihilation. Since in standard co-annihilation, 
all heavier states in the dark matter sector decay to the lightest, pdm = '"^lightest '^j ^-nd 
Eq. ( |3.6D determines the dark matter relic density directly. This equation is identical to the 
Lee- Weinberg case and thus can be solved by an identical technique. 

When semi-annihilations are present, there is simply no way to guarantee the condition 
in Eq. (p.5|). In fact such a condition is usually not satisfied as we will see in the next section, 
and therefore we have to solve the system numerically. 



4. The bbx model 

While the meson-baryon system in Sec. |3.1| is representative of dark matter models exhibit- 



ing semi-annihilation, numerically solving the 0{N'j) coupled Boltzmann equations is not 
particularly enlightening. Here, we study in detail a minimal multi-component dark matter 
model having semi-annihilations among its allowed reactions. This can be obtained from the 
system in Sec. |3.1| by taking Nf = 2 and then "orbifold" identifying the flavor and baryon 
symmetries. 

In this toy system, dark matter is composed of two stable components, a complex scalar 
field X with mass and a vector-like fermion b and b'^ with mass mf,: 

CbbxJree = l^^xl^ + iba^'d^b + iWa^B^b'' - va\^x - [mbb% + h.c] . (4.1) 

We impose a U{1) global symmetry under which 6, b'^, and x have charges +1, —1, and —2, 
respectively. The U{1) invariant interacting Lagrangian for this model is 

jC-bbx,int = [Ki^bW + K2xbb + K2X^b%^ + h.c] + K3X^X0 + «4X^X<^^ - y{(t)), (4.2) 

where we also introduce a real scalar cj) field, the same portal field as in Sec. which we assume 
to be in thermal equilibrium in our calculation and decays to SM states. The potential V{(j)) 
is of the same form as Eq. (|2.2| ) . 

Fermion number conservation always guarantees that b is stable, whereas stability for x 
requires the triangle inequality rule m-^ < 2mb (isosceles triangle). In terms of symmetries, 
the stability of b is guaranteed by a Z2 subgroup of the original U{1) symmetry under which b 
and b'^ goes to minus themselves with all other fields untouched. Likewise, given the triangle 
inequality, the decay of x to other particles is forbidden by a subgroup of the original 
U{1) symmetry, under which x goes to minus itself (and b and b^ pick up a factor of i and 
—i, respectively). The model and its symmetries are summarized in Table ^. 



- 12 - 



Fields 


Spin 


U{\) charge 


Z4 


Z2 


b 


Weyl left 


+1 


+i 


-1 




Weyl left 


-1 


—i 


-1 


X 


complex scalar 


-2 


-1 


+1 





real scalar 












Table 4: Field content and symmetries of a minimal multi-component model with semi-annihilation. 



4.1 Boltzmann equations 

We can now write down and solve numerically the Boltzmann equations for the relic abun- 
dance of h and x and see how semi-annihilations affect the dark matter relic density. In order 
to write down the system of coupled Boltzmann equations, we have to identify the processes 
which change the number of h and x fields. We have to be careful here, since both fields 
are complex and we have an extra degree of freedom (the antiparticle) . We will follow the 
evolution of the number density of the particle only for both species, and once we have their 
relic abundance, we have just to multiply the result by two to take into account also the 
contribution from the antiparticles. We write the Boltzmann equations for the variables 
and n^, and make the assumption rn, = and for a standard WIMP. 

The reactions which change the number of h particles are 

bb^cfxt), bb^xx, bb ^ ^x, H^bx, bx^H, (4.3) 
whereas the ones which change the number of x particles are 



XX 



XX bb, x4> bb, xb ^ (pb. 



(4.4) 



If we assume CP invariance, the collision operators for b(f) 6% and 6% — > bcj) cancel out 
in the Boltzmann equation for b, since they do not change the total number of 6 + 6. The 
Boltzmann equations take the form 

^+3Hn^ = Y,Ci, (4.5) 



where the sums run over all the collision terms. The collision operators for the reactions 
which change the number of b particles are 



C. 



bb- 



^bb-^XX ~ ~ (^'"^e\)bb-^xx 



2 cq2 



n 



eq2 



rih 



en '^x 

nx 



(4.6) 



n 



eq2 



eq 2 '"X 

nx 
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whereas the ones for the reactions which change the number of x particles are^ 

^Xb~^<j)b = ~ {^'^'ccl) xb-^<j)b [^b^x — "-6^^^^] , 
^XX^bb = - {(^Vj.el)xx^bb 



2 2 

Ti^. kTIi. 

X eq2 o 

^b 



(4.7) 



n 



cq2 

eq-n-x — nl 



"6 „ „2 



Once we fix the dark matter masses and limit our consideration to s-wave annihilation, 
there are four free parameters, which we choose to be the following s-wave matrix elements 

^b-b-^^^ = oi, M^^^^^ = /3, Mf,i^^^ = K, M^f,^fh = e. (4.8) 

The diagrams for the reactions in Eq. (^]^) are shown in Fig. ^, and the matrix elements for the 
other reactions can be obtained by crossing symmetry. The first two amplitudes come from 
the standard annihilation to light particles in the final state, the third from the process where 
two dark matter particles of one species are entirely converted to two dark matter particles 
of the other species. This process is phase space suppressed for typical kinetic energies at 
freeze-out. The last corresponds to the semi-annihilation process, which is the new effect we 
take into account. In particular we are interested to see how the relic abundance changes with 
the parameter e. The explicit expression for the thermal averages relevant for the collision 
operators are given in App. 

In order to find the relic abundance of b and x particles, we have to solve a system of 
two coupled Boltzmann equations. As discussed in Sec. |3.3| , it is not possible to sum the 
equations and get a single Boltzmann equation, as one can do in the case of co-annihilations. 
The reason is that the reaction b(j) — t- x<P does not take place at all, since it violates the U{1) 
symmetry of the model as well as fermion number, thus we cannot make any assumptions on 
the relative number density of the two species. We have to solve the system numerically, and 
to do this it is convenient to write this system in dimensionless variables. We define 

Yi = —, X= Xab^cd{x) = ^ %:{<yv^e\)ab^cd{x), (4.9) 

s T H[x = 1) 

where the time variable x is defined with respect to the mass of x- The Boltzmann equations 



^The collision term for x4' ~^ is written by using the opposite reaction, see Eq. (A.£) in App. ^ The 



factor of 1/2 comes from the thermal average over two identical particle in the initial state, which is not 
compensated by any factor of 2 since in this process we get rid of one single x particle. 
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Figure 5: Reactions in the bbx model. All the other ones can be obtained by crossing symmetry. 
The dark matter particles are drawn as double lines, and the coefficients of the s-wave amplitudes are 



given in Eq. (4.8) 



in the new variables are 



dYb _ 


1 


dx 




dY^ _ 


1 


dx 


x'^ 




+ 



V2 ^^cq2 



bb^XX 



Veq2 

_ lb V2 

^b ^eq2^X 



+ A, 



bb^4>X 



-^eq2 
^b Y 

^x 



bb^X(l> 



\^X ^T"^] + \x^bb 

+ A 



X Y' 

X ^eq2-'b 



T^eq2 

ib V _ V2 

ycq -fx ^b 



Xb^4>b 



[YbYx - YbY^^] 



(4.10) 



where the functions Xab^cd also depend on the time variable x. The functions Y^^ are defined 
analogously to the non equilibrium ones in Eq. ( [4.9| ), and they take the form 



9x 45 ^2 



ycq _ 

^ s Qtfs 47r^ 



n 



eq 



9b 45 __2 2 



yeq ^ 

^ s g.s 47r4 



x'K2[x\ 
r'^x^ K2[rx] 



(4.11) 



where we define r = rriiy/m^ and i^2[a^] is the modified Bessel function. 



- 15 - 



K = 0, e = 0, nh^ =^0.1 K = 2, e = 0, nh^ 2:^0.1 




X 

(c) 

Figure 6: Y vs. x evolution in the bbx model for the b particles (blue lines) and the x particles (green 
lines). The masses are taken to be = 0.8 TeV and mty = 1 TeV. We plot the equilibrium values 
(dashed lines) and the numerical solutions (solid lines). The pure annihilation amplitudes are fixed 
to be a = 1.6 and (3 = 0.8, the values of k, e and JIdm^-'^ are shown in each plot. We see that both 
K (6 o X conversion) and e (semi-annihilation) have an effect on the relative and total dark matter 
densities. 

4.2 Numerical results 

In our numerical study we fix the dark matter masses to be niy. = 0.8 TeV and = 1 TeV 
for concreteness, and we vary the matrix elements amplitudes. We consider three main cases. 
For each one, we fix the values of a and f3 (annihilation amplitudes) and study how the 
relic density is affected once we turn on the parameters k (6 -f-)- x conversion) or e (semi- 
annihilation) . 

Case I: a ~ /3. We start from the case where our choice of a and /3 produces an equal 
amount of b and x today, namely a = 1.6 and /3 = 0.8, as shown in Fig. ^a). If we 
turn on the reaction bb — )• XX) ^ ™ Fig. ^(b), the effect is neglegible and the net result 
is a small dilution of b particles, as expected, since b is heavier than x- If we turn 
on semi-annihilations, we get fewer relic particles, and for e = 2 the relic density is 
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a: = 0, e = 0, n /i^ =i 100 /c = 2, e = 0, 11 /i^ 0.1 




X 

(c) 

Figure 7: Same as Fig. ||, now fixing the pure annihilation amphtudes at a = 0.03 and /3 = 1.5. 
Without species changing interactions, such a scenario would lead to over production of b particles, 
but both K (6 O X conversion) and e (semi-annihilation) can be used to achieve the proper relic 
density. 

reduced by a factor of 2 compared to e = 0, as shown in Fig. ^(c). Thus we see that 
semi-annihilations do affect the dark matter relic density. 

Case II: a <C /3. We now consider regions of parameter space where the pure annihilation 
contributions would create very few particles of one kind and many of the other. In 
the absence of the amplitudes e and k, such a parameter region would be excluded by 
overclosure, but the contribution of these species changing reactions can bring the relic 
abundance back to the observed value. When a <C /? the pure annihilations would give 
many b and very few sls shown in Fig. |^(a) for a = 0.03 and /3 = 1.5. If we turn k 
on, as in Fig. ^(b), such a reaction converts b particles into x particles very effectively, 
and the latter have a huge annihilation cross section. The net result is a reduction of 
the total dark matter particles b and Xi and for k = 2 we get the relic abundance we 
observe today. Analogously, when we turn on the semi-annihilation bb — )• (j)x we can 
get the observed relic density, as in Fig. 0(c) for e = 1.5, since this reaction destroys b 
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a: = 0, e = 0, n /i^ =i 100 k = 4 n, e = 0, n =^ 100 




X 

(c) 

Figure 8: Same as Fig. now fixing the pure annihilation amplitudes at a = 2 and /3 = 0.01. Here, 
X particles would be over produced, but because of phase space suppression k (6 O x conversion) is 
ineffective. However, e (semi-annihilation) is never phase space suppressed and can restore the desired 
relic density. 

particles to produce x particles. 

Case III: a ^ (3. From the previous two cases, one might erroneously conclude that 6 -f-T- x 
conversion and semi-annihilation have similar effects on the relic density. However this 
is not the case, since b -H- x conversion is phase space suppressed while semi-annihilation 
is never phase space suppressed. This is beautifully illustrated in the case when pure 
annihilations would give a lot of x arid very few b, as shown Fig. |^(a) for a = 2 and 
/3 = 0.01. The semi-annihilation process xb — ^ b<j3 is still able to destroy x particles 
and bring the relic density back to the observed value. Fig. ^(c) for e = 5. However, 
the reaction xx ~^ bb is now powerless, even if we badly break perturbation theory for 
K = 47r as in plot Fig. |8|(b). This makes sense, since such a reaction is forbidden at 
zero kinetic energy and at freeze-out, the thermal energy is very small compared with 
the mass splitting between the two particles. Thus, semi-annihilation is a truly unique 
species changing interaction that affects early universe cosmology. 



- 18 - 



We have considered only nib > above discussion. The results for rriy. > rrib are 

similar, the main difference being that the conclusions of case II and III are reversed. 



5. Implications for indirect detection 



We have seen how semi-annihilation can affect the predictions for dark matter relic abundance. 
Now we explore the implications of semi-annihilations for dark matter detection experiments. 
Since semi-annihilation does not give any new contributions to direct detection rates, we focus 
on indirect detection only.^ 

There has been much recent interest on indirect detection of dark matter, motivated 
by observations suggesting the presence of a new primary source of galactic electrons and 
positrons. PAMELA reported an unexpected rise with energy of the e^/(e^ + e~) fraction in 
comparison with the estimated background [^] , supporting what already observed by HEAT 
[ ^ ] and AMS-01 |41]. However no antiproton excess was found |42]. Such an anomalous 
excess of the total e+ + e~ flux was confirmed by FERMI Q and HESS [||, ||], although 
the observed spectrum is softer than the peak reported by ATIC |46| and PPB-BETS [47j. 
These anomalies might be evidence for dark matter annihilation or decay in our galaxy |45, 



50| , p\\ , |5^ , though astrophysical explanations have been put forward [gS], U, |55|, q6|, p7|| . 



In general, semi-annihilation does not lead to the "boost factor" necessary to explain the 
magnitude of these candidate dark matter indirect signals. However the presence of such a 
reaction can still have considerable implications for indirect searches, as we shall see explicitly 
for the models considered in Sec. ^ and Sec. ^ Semi-annihilations are additional channels to 
produce light particles in dark matter interactions in our galaxy, thus the predicted spectrum 
is enriched with respect to the standard scenarios. In addition, our examples models can 
account for the fact that dark matter (semi-)annihilates preferably into leptons; if we assume 
a portal (p field mass < 2GeV the decay of (j) to antiprotons is kinematically forbidden 



5.1 Cosmic rays via a scalar portal 

In the class of models we consider, semi-annihilations are new reactions to produce light 
particles in the final state, and for which there is no analog in standard multi-component 
scenarios.^ We now study the injection spectrum of (j) particles in the models we considered 
in this paper. As already said, we focus on semi-annihilation to a portal field <j), but in more 
general dark matter scenarios (p could be a SM field itself. 

'^For an extensive collider study of models with a Z3 symmetry see Ref. |^]. For dual-component WIMP 
dark matter without semi-annihilation in direct detection, indirect detection, and collider experiments see 
Ref. [|||. 

*There are DM scenarios where DM annihilates into two light states, one unstable and one stable jssl ]. 
Strictly speaking, this is an example of semi-annihilation, but because the stable final product is much lighter 
than the DM particles, the kinematics is effectively equivalent to ordinary annihilation. 
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The differential energy flux of 4> from dark matter annihilation and semi-annihilation in 
the Milky Way is of the form 

^ oc n,n, {avU, 6 (e^ - 4) , (5.1) 

^ ij 

where the sum runs over the dark matter particles, and Nf^ and E^^ are the number and the 
energy of (j) produced in any reaction, respectively. Thus, the (p spectrum consists of a number 
of monocromatic lines with different intensities. The spectrum of the produced SM particles 
which we observe in indirect detection experiments would be obtained by convoluting the (f) 
spectrum with the interactions in Eq. ( [2.21 ). We consider here as an example the particular 
case where (p has a two body decay (p — ^ 77- The decay products are isotropic in the (j) rest 
frame, therefore the energy distribution in the galactic frame is flat. If we assume the mass 
hierarchy moM ^ n^(i> ^ "t-^i the energy distribution of 7 in the galactic frame has support 
< < E^^ for any reaction, and its differential flux takes the form 

5.2 Spectra in Z3 model 

We start with the model with Z3 symmetry considered in Sec. ^, where we have^ 

oc nl [{av)^^^^^ + {av)^^^^^] . (5.3) 

First, we note that the integrated flux of (p particles does not differ from the standard case 
where the thermal production is dominated by annihilations only, namely oc nQ{av)o. 
The relic density constraint imposes n^. = uq, and it also imposes the equality between the 
{crv)o and {av)^^^,j)^ + (av)^^^^^ (see Eq. (^!7|)), therefore = '^'^a- Although semi- 
annihilations cannot provide us with any boost in the indirect detection rate, they do change 
the injection spectrum of (p particles. 

In the pure annihilation case, all the (p particles from dark matter annihilation are pro- 
duced with energy equal to the mass rriy. of the x particle, but in the semi-annihilation 
reaction, (p in the final state has energy 3m^/4. In the most general case we have both contri- 
butions, thus the (p injection spectrum has two monochromatic lines at 3m^/4 and m^, and 
the relative amplitude of such peaks are given by the ratio of the squared amplitudes. These 
two quantities are related by the relic density constraint found in Sec. 4?7^ -|- 3r]^ ~ 1.5, 
thus, to a good approximation, the two ratios depend on only one amplitude, namely 

where is the flux of (p particles with energy E. The results for the (p spectrum are 
shown in Fig. In the limiting case where only one process is present, the (p spectrum is 

^The semi-annihilation produces only one particles, however there are two contributions, from XX ^-nd 
XX ill the initial state, thus there is no relative factor of 1/2. 



- 20 - 



]]a = 0.6, T]s = T]a = 0, l]s = 0.7 




E/m ns 



(c) (d) 

Figure 9: Injection spectrum of (j) particles in the model with symmetry for: (a) annihilation 
only; (b) semi-annihilation only; (c) equal amplitudes from the two processes. The total spectrum in 
normalized to 1. The relative line intensities as a function of rjg are shown in (d), i?^ (blue line) and 
^tm/i (yoUow line) . The results are independent of the particle mass m-^ . 



monochromatic, see Fie;. 1(a) and Fig. |(b). When both processes contribute, the spectrum is 
in general bipolar, as shown in Fig. |9|(c) for the case of equal amplitudes for the two reactions. 
The relative spectrum amplitudes in Eq. ( ^.41) as a function of r/^ are shown in Fig. |9|(d). In 
the case of only one relic particle, the results are independent of the mass (up to a small 



logarithmic correction through the freeze-out value in Eq. (2.6)). 
5.3 Spectra in bbx model 

We now analyze indirect detection in the bbx model, which exhibits a much richer kinematic 
structure. In this case, there are four reactions to produce (p particles in the final state, the 
ordinary annihilations of the two components and two semi-annihilations. The overall 4> flux 
now takes the form 

<^bb^ oc nl [{crv,ci)bb-,^^ + {crVrci)bb^cPx] + ^brixicrv^^l) bx^b^p + nl{av,el)xx^<l><P- (5-5) 

Let us first briefly discuss the possibility of achieving a boost factor in this dual-component 
model. We numerically evaluated the flux in Eq. (|5.5|) for the region of parameter space which 
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gives the correct relic abundance, and compared such a flux with the one obtained when only 
the lightest particle is present (this is a conservative comparison, since the flux in the latter 
case is proportional to m^^). We found that the most we can get is a boost factor of a 
few in the total indirect detection rate. Of course, semi-annihilation could coexist with an 
alternative boost mechanism such as the Sommerfeld enhancement ISOl, |60|, IsTtl or a local dark 



matter overdensity pU 62, 63 



The spectrum is this multi-component case is quite rich, since there are now four processes 
giving (p in the final state, yelding four monochromatic (p lines. The energies of the four lines 
are given by the energies of (j) in the final state 

rn^ 2mi, + m-^ 

Ef,i^^^ = mb, Ebb^^^ = nib - ^bx^H ^ ^ rub + ' ^xx^<P4> = ^x- (^-^J 

We present numerical results for the (j) spectrum for the same values of the dark matter 
masses chosen in Sec. namely m^^ = 0.8 TeV and nib = 1 TeV. This gives E^^i^^^ = 1 TeV, 
Ebb^<f>x = 0-84 TeV, E^^^i^ = 0.62 TeV and E^^^(j,(j, = 0.8 TeV. To simplify the following 
discussion we perform our analysis for k = 0. The inclusion of a finite value of k does 
not change the expression in Eq. (5^), since its associated process does not produce any cj) 
in the final state, but it can change the thermal relic density as discussed in Sec. ^ Two 
example spectra are shown in Fig. 10, where two opposite cases are considered. In the first 
case, Fig. |To|(a), the semi-annihilation contribution is taken to be very small, thus the pure 
annihilation lines E^^i^^^ and E^^^^j,^ dominate. When the thermal production is mostly 
controlled by semi-annihilations, as in Fig. |l^(b), the semi-annihilation lines Ebb^4,x ^"^^ 
^hx-^H overwhelm the standard contributions. 

To study more carefully the relative intensities of the four lines, we fix the ratio q//3 = 2 
(as in the two spectra just shown) and keep k = 0, leaving us with two free parameters, a and 
e. We take the relic abundance constraint, J^dm^^ — 0.1, to get a as a function of e, and we 
study how the line intensities change as we increase the contribution of semi-annihilations. 
The result is shown in Fig. p!0|(c). For very small e, the lines from semi-annihilation give a 
negligible contribution, whereas the standard lines equally share the overall fiux (since for 
a = 2/3 and small e we get approximately as many 6 as x)- Once e gets big enough (e ~ 1) the 
thermal production is controlled as much by annihilations as by semi-annihilations, therefore 
the four lines have approximately equal intensity. For bigger values of e the annihilation 
lines intensities are very suppressed, and eventually for e > 3 the line with energy E^^_^i^ 
dominates the signal. The reason why for large e only one line is present is easily understood 
looking the Boltzmann equations system in Eq. ( [4.10 ), and considering only the collision 
terms associated to semi-annihilations. The collision operator associated to the reaction 
X(t> — ^ bb in the equation for Yy- tends to increase the total number of thus the dark matter 
density balance today is totally dominated by the x particles, and the line associated to the 
semi-annhilation with x iii the initial state gives the dominant contribution. 

One expects the indirect detection spectra in the meson-baryon system in Sec. 3J to have 
an even richer structure of (j) lines. Of course, once convolved with the (f) decays, the spectrum 
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Figure 10: Spectra of particles in ttie hhx model for two opposite cases, namely thermal production 
mostly controlled by: (a) annihilations, (b) semi-annihilations. The values of the matrix element 
amplitudes are shown in each plot. The relative intensities of the four lines are shown in (c) as a 
function of e for a — 2(3 keeping the dark matter relic density fixed: Ej^^^^^ (blue solid line), Etb^tpx 
(red dashed line), E^^^^i^ (green dashed line) and E^^^^^ (yellow solid line). We always take k = 0, 
= 0.8 TeV, and mi, = 1 TeV. 




is less distinct. The injection spectrum when — )■ 77 (see Eq. ( p.2[ )) is shown in Fig. for the 
same two cases considered in Fig. |l^(a) and Fig. |l^(b). As we increase the semi-annihilation 
contribution, the differential energy spectrum gets shifted to lower energies. This result is 
consistent with Fig. 10(c), since bx — ^ bcj) dominates for large e, and such a process has the 
less energetic line. More detailed study is necessary to know whether such a spectrum can be 
measured in 7 rays telescopes. 
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Figure 11: Differential plioton energy spectrum in the bbx model for the case — 77. The same 
two cases as in Fig. ^are considered, and the overall signal is normalized to 1. 



6. Conclusions 

The presence of dark matter in our universe motivates us to look for physics beyond the 
SM, and understanding its origin and composition is one of the great open questions in 
particle physics and cosmology. Particle physics models explaining the origin of the Fermi 
scale often contain stable massive particles called WIMPs. Assuming WIMPs are thermally 
produced in the early universe, these particles have the right annihilation cross section to give 
the observed relic abundance, Odm^^ — 0.1, making the WIMP paradigm one of the best 
solutions to account for dark matter. 

In this work, we relaxed some assumptions about the symmetry structure of WIMP 
interactions by allowing relic particles to "semi-annihilate" . We saw that the semi-annihilation 
process 

Vfc0 (6.1) 

can have a considerable effect on the dark matter relic abundance. Such processes are present 
when dark matter has a larger stabilizing symmetry than just Z2, and we showed explicit 
examples where semi-annihilation is present: a single species dark matter model with a Z3 
symmetry and a multiple species dark matter model with "baryon" and "flavor" symmetries. 
Along with the standard WIMP scenarios summarized in Table 0, semi-annihilation does 
occur in realistic particle physics scenarios. Indeed, as sketched in App. the simple case of 
a supersymmetric QCD theory with Nf = + 1 generically involves semi-annihilation. 

We saw that when dark matter is composed of just one species, the effect of semi- 
annihilation on the relic density can be derived through a standard freeze-out calculation 
[^. Semi-annihilation can in fact dominate over ordinary annihilation for some regions of 
parameter space. However, such single species models are not representation of generic semi- 
annihilating dark matter models, which usually involve more than one stable dark matter 
component. In such cases, a semi-analytical solution for the relic density is simply not pos- 
sible, and one must resort to numerically solving the complete set of coupled Boltzmann 
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equations given in Eq. (|3.3| ). In particular, the simplifying assumptions made when analyz- 
ing co-annihilating models |jl^ are not applicable for semi-annihilation. As a consequence, 
the dark matter dynamics in the presence of semi-annihilation are far more varied than in 
decoupled multi-component scenarios. 

While the inclusion of semi-annihilation does not give any new contributions to the 
direct detection of dark matter, it does have interesting implications for indirect detection 
experiments. We studied the injection spectrum of light particles in the two toy models 
we studied, where dark matter (semi-)annihilates through a (j) "portal" which subsequently 
decays to SM states. The overall integrated flux from (semi-)annihilates in the Milky Way 
halo is not very affected by semi-annihilation (in the case of one relic particle, is not affected 
at all). However, the final state spectrum in semi-annihilating models is far richer than 
in standards scenarios because of the differing kinematics between semi-annihilation and 
ordinary annihilation. 

In the language of Ref. |15|, semi-annihilation is in some ways the "fourth exception" in 
the calculation of dark matter thermal relic abundances. We find it intriguing that unlike 
the traditional three exceptions (co-annihilation, annihilation below a mass threshold, and 
annihilation near a pole), this fourth exception not only affects dark matter interactions in 
the early universe, but also leaves an imprint today via the indirect detection spectrum. We 
expect that there are a wide variety of multi-component dark matter models with species 
changing interactions, motivating further studies of the semi-annihilation process. 
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A. Boltzmann equations in FRW background 

In this appendix, we briefly review the derivation of the Boltzmann equations in an expanding 
universe. If the relic particle, which we assume to be massive, remained in thermal equilib- 
rium until today, its abundance would be negligible because of the exponential factor in the 
Maxwell-Boltzmann distribution. But if the interactions which keep the particle in thermal 
equilibrium freeze-out when the temperature is not small compared with the relic particle 
mass, they can have a significant density today. 

We consider a generic particle ij^a whose number density evolution is governed by the 
Boltzmann equation 

^ + 3Hna = Ca, (A.l) 

where H is the Hubble parameter and Ca is the collision operator for the reactions which 
change the number of ipa particles. The 3Hn factor accounts for the dilution effect of the 
expansion of the universe, and the right hand side accounts for the interactions that change 
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the number of ip. If there are no interactions, the right hand side is zero and we have 
UaR^ = const, where R is the scale factor, simply due to the expansion of the universe. 
Considering only 2^2 processes, we can write the collision operator as 

Ca = ^ Cabled, (A. 2) 

b,c,d 

where the sum runs over all the possible reaction ab — > cd. The collision operator for a single 
reaction results in 

Cab-^cd = - J (27r)^(5^ {Pa+Pb -Pc- Pd) dRadJlbdncdnd \Mab^cd? fafb - IMcd-^abl^ fcfd 

(A.3) 

where fi is the phase space distribution for the particle i. The relativistic invariant phase 
space is defined as 

where Qi is the number of the internal degrees of freedom of the particle i. If we assume CP 
(or T) invariance, we have an equality between the two matrix elements Mab-^cd = •^cd-^abi 
and we can rewrite 

Cabled = - j (27r)^5^ {pa + Pb-Pc- Pd) dYladnbdncdUd \Mab^cd\'^ [fafb - fcfd] ■ (A.5) 
In the following, we will also assume kinetic equilibrium during freeze-out 

f,iE,t) = ^fP{E,t), (A.6) 
with the equilibrium number density being 



-r=5. /7^r(p), (A.7) 



d^p 
(2^ 

where /f^(p) is the thermal distribution for the species when it is in thermal equilibrium. 
Thus the spectrum remains thermal up to a constant which takes into account the fact that 
the relic particles are disappearing. The thermal distributions satisfy the following relation 
(by energy conservation) 

/aVr = exp [- (E, + Eb) /T] = exp [- (E, + E,) /T] = f^Vf. (A.8) 
By using these two results and performing the integrals over the momenta, we get 



Cabled — (•^'I'reOab— >cd 



en en 
n-alT-b en ea^cf^d 



-{av) 



cd-^ab 



eq eq 

eq eq^a^fe '^c'f^d 



(A.9) 



where a is the total cross section for that particular process, Vj-ei is the relative velocity, and 
{(TVrei) ab^cd denotes the thermal average 

/ \ _ / d^Pad^Pb a (Pa, Pb)ab^cd ^rcl fa'^{Pa)fb'^{Pb) 

!d^Pad^Pbfl\pa)fl\Pb) • ^^-^^^ 



- 26 - 



In Eq. ( A. 9 ) we have two equivalent expressions for the colhsion operator associated to 
the reaction ab — )• cd. The first one is the usual form involving the thermal averaged cross 
section for that particular reaction, while the second one contains the thermal averaged cross 
section for the opposite reaction, cd — t- ab. We can use either form to write the collision 
operators, and we use this freedom in the bbx model, since it is convenient to choose the 
reaction in which the massless degrees of freedom are in the final state in order to have a 
simpler expression for the thermal average. 



B. Supersymmetric QCD example 

In this appendix, we study a toy dark matter model based on supersymmetric QCD where 
the number of flavors A'^^ and the number of colors are related by Nf = Nc + 1 (for a model 
of supersymmetric QCD dark matter where Nf < see Ref. [^). In the ultraviolet, the 
degrees of freedom are quark and antiquark flelds, and Qj respectively. To have a viable 
phenomenological model we also introduce a mass term in the superpotential 

WvY = mlQ'Q-. (B.l) 

In the absence of this superpotential, the theory has a global U{Nf)L x U{Nf)fi symmetry, 
which we use to put the mass matrix in a diagonal form, ml = rriidl. 

Below the QCD dynamical scale, the appropriate degrees of freedom are the meson and 
baryon color singlet moduli, interacting via the superpotential 

WiR = BjMuBi - c^f A^-^f det (M) + c ATr (mM)l , (B.2) 

167r^ 

where we invoked "loop-democracy" and the c factor comes from canonically normalizing the 
meson fields. Using naive dimensional analysis |67] we have c ~ Att. The supersymmetry 



preserving minimum is found by solving the system of equations 



dBi 



^ c^B.M,;,- = 0, 



dWiK 



= c3 MjiBi = 0, 



[c^ BjBi - c^f A^-^f cofMji + c A niij] = 0, (B.3) 



where the vacuum expectation value (vev) of a generic field $ is denoted by boldface type 
The cofactor matrix is defined as cofMji = {—y~^^Aji, where Aji is the minor of the matrix 



Mji. 



The first two equations can be seen as homogeneous linear system with the baryon and 
antibaryons fields as unknowns and the meson matrix vev as the coefficient matrix. We 
assume that det (M) 0, and then the systems in the first two lines of Eq. (|B.3| ) have only 
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the trivial solution B,- 



Bi 



0. Since the matrix M is invertible we can further simplify 



the vacuum condition for the meson matrix. The inverse of M is given by 



M 



-1 



1 



(cofM)' 



1 



det(M) 

and in conclusion the vacuum of the theory results in 

1 



det(M) 



cofMjj, 



B = B = 0, 



M 



=^ n 



mi 



,4=1 



diag I , m. 



1 

2 ' ' 



(B.4) 



(B.5) 



The meson matrix vacuum expectation value is diagonal, thus we have an unbroken 
accidental C/(l)^^ flavor symmetry, as well as a U{1)b symmetry acting on the baryons. The 



spectrum of the model is easily obtained by expanding the superpotential in Eq. ( B.2 ) around 



its minimum, given in Eq. (B^), to quadratic order: 



det(m) 



Nf-l 



BiBi + 



mi 



det(m) 



1 

JVf-l 



E 



mi mi 



A 



(MijMji - MiiM. 



(B.6) 



The baryon and the off-diagonal mesons fields are already mass eigenstates, whereas the 
diagonal mesons have mixing terms. In the limiting case where the mass matrix in the UV 



superpotential in Eq. (BJ.) is a multiple of the identity with all the masses equal to m, we 
have the following scaling for the masses of the IR degrees of freedom 

JVf-2 -I -, JVf-2 

(B.7) 



rriB 



rriM 



An example spectrum is shown in Fig. 12 for the case Nf = 4 and for same particular values 
of the mass mj. In general, supersymmetry breaking effects will lead to splitting within the 
chiral multiplets. 

As already said, there in an unbroken U{l)^f xU{1)b symmetry which makes the lightest 
baryon and the lightest off-diagonal meson absolutely stable, thus they are good dark matter 
candidates. The heavier baryons and off-diagonal mesons might also be stable if their decay is 
kinematically forbidden, as already discussed. In general a few of them can be unstable, but 
it is quite common to have more than just two stable components, and in the limiting case 
when the mass matrix m is almost degenerate, all of them are stable. The diagonal mesons do 
not carry any flavor quantum numbers. In analogy with the meson-baryon system in Eq. ( p.lD 
we denote them by Ma = (pi, and their possible interactions with SM fields would come from 



the supersymmetric analogs of Eq. (2^). From expanding the superpotential in Eq. 
one sees that semi-annihilations reactions involving the stable dark matter particles can take 
place at freeze-out, as for example 



M12M13 -)■ M23</)1, B2M12 Bi(j)i 



(B.8) 



where M, B, and (j) can refer to either the fermions or scalars within the respective chiral 
multiplets. Therefore the standard treatment for the relic density calculation cannot be 
applied. 
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Figure 12: Supersymmetric QCD spectrum for TV/ = A^c + 1 = 4 and for particular values of the 
input masses mi. 

C. Thermal averages in the hhx model 



In this appendix, we derive expressions for the thermal averaged cross sections from Eq. ( A.IOD 
for the hbx model. We work in the s-wave limit where the cross section does not depend on 
external momenta, using the amplitudes defined in Eq. (4 



If the process ab — t- cd is not kinematically suppressed, then in the s-wave limit, the 
integral over momenta in Eq. ( A.10| ) is trivial and simplifies to 



{av)ab^cd = [o- (Pa, Vb)ab-^cd ^rcl]p^=p^=o 



(C.l) 



It is convenient to define 



Wab~,cd{s) = EaEb (Jab-^cdVrel = \ j \^ab^cd\^ d^^'^^ {s,mc,md) , (C.2) 

where s = {pa and d^^'^^ {s,mc,md) is the relativistic invariant two-body phase space. 

In the collision center of mass frame, we have 



s,mc,md) 



dn 

32^ 



(C.3) 



where d^ is the differential solid angle element. Plugging this expression into Eq. ( p.2[) , and 
assuming that the matrix element has no angular dependence (as is the case in the s-wave 
limit), we find 



Wab^cd{.s) 



\M.ab^aA 

327r 



~^ 2 



(C.4) 
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Thus, the thermahy averaged cross section is 



{ov) 



ab—>cd 



En Eh 



Wab^cd (("T-a + mfe)2) 



Pa=P6 = 



rriamh 



(C.5) 



We can use Eq. ( p.5| ) to compute the thermal averages for ah the reactions in the bbx 
model, except for the reactions bb -f-T- xx where the presence of a kinematic threshold requires 
more care |15|. For future reference we also define r = mh/m^, and we give all the results in 
terms of rriy. and r. We have three different reaction types: annihilation, semi-annihilation, 
and dark matter conversions. 



Annihilations 

The ordinary annihilation reactions are bb 
effectively massless final states, thus 



and XX 



1 



a 



A 



i. In this case, we have two 



(C.6) 



where the factor of 1/4 comes from the average over the spin degrees of freedom, and 

/32 



{av 



rel/xx- 



327rm| 



(C.7) 



Semi-annihilations 



In the semi-annihilation case, we have only one massless particle in the final state, thus there 
is an additional phase space suppression 



1 



47-2 



1 



4r2 / 327rm^ 



and 



(o-frel) 



Xb—^<l>b 



2r 



(C.8) 



(C.9) 



Dark matter conversions 



In this last case, we have two massive particles in the final state. The phase space now is 
suppressed, and the reaction can be forbidden at zero relative velocity. To correctly compute 
the thermal average we consider the second exception discussed in Ref . [ p!5| ] , namely when the 
reaction takes place near a kinematic threshold. 

We start from the reaction bb — )• xX- The Mandelstam variable in the center of mass 
frame of the collision is 



s = {E^ + E^f = 4mlj^ 



Ami 



1/2 



1/2 



V2, 



(C.IO) 
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where V2 is the velocity of the x particles in the final state in the center of mass frame for 
the collision. In the s-wave limit we can write 



(C.ll) 



where the factor of 1/4 is because the average over the spin degrees of freedom. We define 

/ib = |r2-l|^/^ (C.12) 



and for the thermal average we have the following expression [15 



1 2 ,,2 1/2 



b^xx 4^2 rV2vrV 

A. 



(C.13) 



where the + sign (— sign) is for the case when the annihilation if allowed (forbidden) at zero 
relative velocity in the initial state. 

For the opposite reaction, XX ~^ bb, the discussion is analogous, up to some factor. We 
do not have to average over any initial spin, and we can write 



327r 



(C.14) 



where V2 is the velocity of the b particles in the final state in the center of mass frame for the 
collision. We define 

1 



1 



1/2 



and for the thermal average we apply again the results of Ref. [15| 



327rm| vr^/ 



(C.15) 



(C.16) 



where it for the allowed/forbidden reaction at zero relative velocity in the initial state, re- 
spectively. 
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